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We investigate the dynamics of neutral atoms in a 2D optical lattice which traps two distinct 
internal states of the atoms in different columns. Two Raman lasers are used to coherently transfer 
atoms from one internal state to the other, thereby causing hopping between the different columns. 
By adjusting the laser parameters appropriately we can induce a non vanishing phase of particles 
moving along a closed path on the lattice. This phase is proportional to the enclosed area and we 
thus simulate a magnetic flux through the lattice. This setup is described by a Hamiltonian identical 
to the one for electrons on a lattice subject to a magnetic field and thus allows us to study this 
equivalent situation under very well defined controllable conditions. We consider the limiting case 
of huge magnetic fields - which is not experimentally accessible for electrons in metals - where a 
fractal band structure, the Hofstadter butterfly, characterizes the system. 



I. INTRODUCTION 



The recent experimental progress in manipulating and 
controlling trapped neutral atoms in optical lattices by 
quantum optical means T, 2] allows for a number of novel 
applications in a variety of different fields like quantum 
information processing |3i iJfe iSi iL »Sj j atom interfer- 
ometry Q , and atomic and molecular physics 0, 0, 0| . 
One of the most important features in all of these ap- 
plications is the large degree of control by quantum op- 
tical techniques over the structure and the parameters 
of the Hamiltonian describing the atomic system. This 
control allows to realize and deploy a number of lattice 
Hamiltonians T^, which are frequently used as toy 
models for strongly correlated condensed matter systems 
and therefore theoretically very well studied. However, 
many of the most interesting effects in strongly corre- 
lated 2D systems appear if an external magnetic field is 
applied Q. Apart from rotating an atomic cloud [T5| . 
as assumed for example in the study of Laughlin states 
with bosonic atoms in Ref. there seems to be at 

present no obvious way of implementing lattice Hamilto- 
nians resembling the effects of magnetic fields with neu- 
tral atoms. 

In this paper we propose a 2D setup for neutral 
atoms which allows to engineer terms in single-band Hub- 
bard Hamiltonians corresponding to an external mag- 
netic field. An optical lattice provides the discrete pe- 
riodic spatial structure and we will use lasers instead of 
a magnetic field to induce a phase for particles hopping 
around a closed path in the lattice resembling an effective 
magnetic field. We will show that the strength of this ef- 
fective magnetic field can be varied by laser parameters 
and be made arbitrarily large, a situation first theoret- 
ically investigated by Hofstadter jlJI for electrons. The 
studies by Hofstadter predicted the emergence of frac- 
tal energy bands e that resemble the shape of a butter- 
fly when plotted against the parameter a = ABe/iirTi, 
where A is the area of one of the elementary cells of the 
lattice, B is the strength of the magnetic field, and e is 
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FIG. f: The Hofstadter butterfly: The eigenenergies e are 
shown as black dots for different a, and J is the hopping 
energy as defined in the text. The most dominant feature is 
the splitting of the energy band into r subbands for a = 1/r. 
For further details of this plot see [l^ . 



the charge of the particles (cf. Fig ^ . The phase 2TTa 
is gained by the wave function of a particle due to the 
magnetic field when it hops around a plaquette of the 
lattice. As shown by Hofstadter the nature of the energy 
bands depends crucially on the parameter a. If a = p/r 
with p, r integers, i.e. a is a rational number, the energy 
spectrum splits into a finite number of exactly r bands 
whereas if a is irrational the energy spectrum breaks up 
in infinitely many bands and thus the fractal structure 
shown in Fig. ^ emerges (for a detailed discussion on the 
properties of the Hofstadter butterfly see jU]). 

Fractal energy band structures are believed to play an 
important role for a number of effects like the quantum 
Hall effect induced by mag netic fields in strongly corre- 
lated electron systems [ij, • Therefore it is desirable 
to study the Hofstadter butterfly experimentally under 
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well defined conditions. However, it turns out that the 
area of the elementary cells in metals where fractal energy 
bands could possibly be seen is so small that huge mag- 
netic fields would be required to obtain values of a which 
are on the order of 1 [13 • Also in more sophisticated su- 
perlattice setups with larger area A it is experimentally 
very difficult to obtain direct clear experimental evidence 
of the Hofstadter butterfiy 19]. 

We will consider a 2D system of ultracold atoms 
trapped in one layer in the xy plane of a three dimen- 
sional optical lattice. The atoms are in the lowest mo- 
tional band which can be achieved e.g. by loading the 
optical lattice from a Bose-Einstein condensate [l|, |23| • 
Hopping along the z direction is turned off completely 
by the lattice potential. Different columns of the lat- 
tice trap atoms in internal states! 0) (|e)) (denoted in 
Fig. by open (closed) circles) |2, 01 ■ addition the 
optical lattice is either accelerated along the a;-axis or 
an inhomogeneous static electric field is applied and two 
Raman lasers driving transitions between the states \g) 
and |e) induce hopping along the a;-axis while hopping 
along the y-axis is controlled by the depth of the optical 
lattice along this direction. This setup corresponds to 
applying a magnetic field with a parameter a — qX/An 
where q is the wave number of the Raman lasers along 
the y direction and A the wave length of the lasers cre- 
ating the optical lattice. We also note that the atomic 
setup we are going to describe here can be used for a 
large number of other purposes. It is straightforward to 
add terms to the system that correspond to an external 
electric field. Also the atoms will interact with each other 
via collisional contact interactions 'F| and off site interac- 
tions can be engineered by dipolar Rydberg interactions 
1^. Therefore the setup presented here can be used for 
a number of studies related to the behavior of charged 
particles in a 2D configuration subject to magnetic and 
electric fields and also to study strongly interacting and 
thus strongly correlated systems. Furthermore it might 
be possible to extend this model to different geometries 
of optical lattices. 

In this work we will concentrate on a possible setup 
required to implement the effective magnetic field in an 
optical lattice. We will discuss in detail the laser setup 
which leads to an effective magnetic flux through the op- 
tical lattice, and calculate the corresponding matrix ele- 
ments in Sec.m We also show that it is possible to reach 
each point within the Hofstadter butterfly apart from a 
negligibly small region around a = with the proposed 
setup. In Sec. lIIII we suggest one possibility of measuring 
some of the basic properties of the Hofstadter butterfly 
and discuss the limitations on the resolution for measur- 
ing the energy bands. We also give a brief account of 
the interaction effects. Finally we conclude with a short 
outlook on how the present setup could be extended in 
Scc. lIVI While the focus of the present work is the deriva- 
tion of the single-particle terms in the Hubbard Hamilto- 
nian mimicking a strong magnetic field, we see as one of 
the main motivations the extension of strongly correlated 



many-atom systems in strong (effective) magnetic fields. 



II. SETUP AND MODEL 

In this section we discuss the experimental setup re- 
quired to produce a Hofstadter butterfiy for neutral 
atoms. We first present the optical lattice setup, then in- 
troduce an additional acceleration or static electric field 
and finally describe in detail the additional lasers re- 
quired for our purpose. 



A. Optical lattice 

We consider a three dimensional optical lattice created 
by standing wave laser fields which generate a potential 
for the atomic motion of the form (we use h = 1 through- 
out the paper) 

V{x) = Vox sin^(fca::) + Voy sin^ {ky) + Voz sin^(fcz), (1) 

with k — 27r/A the wave- vector of the light and spatial 
coordinate x = {x,y,z}. The recoil energy is given by 
Er = fc^/2M with M the mass of the atoms. We as- 
sume the lattice to trap atoms in two different internal 
hyperfine states |e) and \g) and the depth of the lattice 
in X- and z-direction to be so large that hopping in these 
directions due to kinetic energy is prohibited |2(]| |. Fur- 
thermore we assume that adjusting the polarization of 
the lasers which confine the particles in the a;-direction 
allows to place the potential wells trapping atoms in the 
different internal states at distances A/4 with respect to 
each other |3, 0, ^ts shown in Fig. QJl. Therefore the 
resulting 2D lattice has a lattice constant (disregarding 
the internal state) in x-direction of a^, = A/4 and in y- 
direction oi Uy — A/2. We restrict our analysis to one 
layer of the optical lattice in the xy plane since in the 
following there will neither be hopping nor interactions 
between different layers. The dynamics of bosonic atoms 
occupying the lowest Bloch band of this optical lattice is 
wefi described by the Bose-Hubbard model (BHM) 2^ 

n,m 

rtGeven,m n,m 

'^Ux ^n,m'^n+l.m^'n+l,mO'n,m: (2) 
n,m 

where is the hopping matrix element for particles 
to tunnel between adjacent sites along the y-direction. 
The energy difference between the two hyperfine states 
is uJeg > and the operators an,m (aJi m) bosonic 
destruction (creation) operators for atoms in the lowest 
motional band located at the site which is centered at 
x„,m = {xn,ym}, whcrc x„ = nA/4 and = mA/2. 
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FIG. 2: Optical lattice setup: Open (closed) circles denote 
atoms in state \g) (|e)). a) Hopping in the y-direction is due to 
kinetic energy and described by the hopping matrix element 
being the same for particles in states |e) and \g). Along 
the a;-direction hopping amplitudes are due to the additional 
lasers, b) Trapping potential in a;-direction. Adjacent sites 
are set off by an energy A because of the acceleration or a 
static inhomogeneous electric held. The laser fii is resonant 
for transitions \g) ^ \e) while 0,2 is resonant for transitions 
between |e) « \g) due to the offset of the lattice sites. Because 
of the spatial dependence of f2i,2 atoms hopping around one 
plaquette get phase shifts of 2Tia — —fm + O + ipm+i+0 where 
ifim ~ mq\/4-K as indicated in a). 



The corresponding mode functions are the localized Wan- 
nier functions w(x — Xn^m) 113 found by suitable super- 
positions of the Bloch functions for the lowest Bloch band 
of the lattice. Since in x-direction the separation be- 
tween two neighboring atoms is half the original lattice 
constant A/2 the overlap of the mode functions of par- 
ticles in adjacent lattice sites might lead to significant 
nearest neighbor interactions described by Ux whereas 
we neglect any other offsite interactions |23|. The pa- 
rameter Un describes the onsite interaction between two 
particles occupying the same site. This onsite interaction 
may depend on the column index n because of different 
internal states with different scattering lengths occupy- 
ing different columns. Since for even [odd] n the op- 
erator an^m describes atoms in internal states |e), [\g)] 
and the Wannier functions for sites which are separated 
by multiples of the original lattice constant A/2 are or- 
thogonal to each other we find the commutation rela- 
tions [an,m,ali^rn'] = ^n,n'Sm,m' with dn,n' the Kroucckcr 
delta. The details of the derivation of the above Hamil- 
tonian can be found in where also the definitions for 
the parameters U, Ux, and are given. 



B. Acceleration or static electric field 

In addition to the above setup we assume an energy 
offset of A between two adjacent sites in x-direction 
as shown in Fig. |2|3. This can be done by accelerat- 
ing the optical lattice along the x-axis with a constant 
acceleration a^cc leading to an additional potential en- 
ergy term for one atom of -ffacc = Ma^ccX- Alterna- 



tively, if both of the internal atomic states |e) and \g) 
have the same static polarizability /i an inhomogeneous 
static electric field of the form E{x) = 5 Ex can be ap- 
plied to the optical lattice yielding a potential energy 
term iJacc — fJ-SEx. We keep this additional poten- 
tial energy small compared to the optical lattice poten- 
tial and treat i?acc as a perturbation. In second quan- 
tization this yields i/acc = ^J2n,m''^^li,m<^n,m whcrc 
A = ^5 EX /A in the case of an inhomogeneous electric 
field and A = Ma^ccX/A if the lattice is accelerated. The 
condition for this perturbative treatment to be valid is 
A <C i^a; with Vx — \/AEjiVox the trapping frequency of 
the optical lattice in the x-direction. 

C. Additional lasers 

Finally, we want to induce hopping along the x- 
direction by two additional lasers driving Raman tran- 
sitions between the states \g) and |e). Each of them con- 
sists of two running plane waves chosen to give space 
dependent Rabi frequencies of the form 

r!i,2 = l^e±*«^ (3) 

with Q the magnitude of the Rabi frequencies, and de- 
tunings ±A. We choose the parameters > 0, A > 0, 
and g > 0. As discussed in detail in Appendix IXI this can 
always be achieved by superimposing two running wave 
laser beams incident in the xy-plane for q > (A + uj^g)/ c 
with c the speed of light. We assume the lasers not to 
excite any higher lying motional Bloch bands and also 
no transitions with detunings of the order of A, i.e., 
i7 ^ A ^ Vx- Then the lasers ^i{2) will only drive 
transitions n — 1 ^ n if n is even (odd) and we can ne- 
glect any influence of the nonrcsonant transitions. We 
find the following Hamiltonian describing the effect of 
the additional lasers 

Hus = (7K,maLm»n-l,m + ^-c) 

rn.n 

-A^nal„^an,m, (4) 

where we have neglected all other terms due to being 
nonrcsonant and defined matrix elements 7„ for even n 

7«,m = ^ y"rf^a; w*(x-x„,„)r2iw(x-x„_i_m), (5) 

and for odd n 

ln,m = ^ j d\ W*(x - X„^„j)f72w(x - X„_i^„J. (6) 

For an optical lattice potential of the form Eq. 

the Wannier functions can be written as a product 

of three normalized Wannier functions, i.e. w(x) = 
'w{x)'w{y)w{z) (23, and we can write 

7n,™ = ^e2"""rj,(a)r,, (7) 
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a) X)/Er b) Vo/Er 

FIG. 3: Matrix elements: a) Fx as a function of the depth 
of the optical lattice Vo/Er. b) Matrix element Fy(Q) as a 
function of the depth of the optical lattice Vb / Er for different 
values of a = {1/8, 1/4,3/8, 1/2} decreasing with increasing 
a. The inset shows Ty{a) against a at Vq/Er = 10. 



atoms in optical lattices and in particular show that it 
can be used to study the whole of the Hofstadter but- 
terfly shown in Fig. ^ We note that wc have chosen the 
detunings of the additional lasers rii_2 to exactly cancel 
the term i?acc arising from the acceleration or electric 
field. If these two terms did not cancel the remaining 
terms resemble a homogeneous electric field. 

III. DISCUSSION 

In this section we show how to detect basic properties 
of the fractal energy spectrum in an experiment. We also 
briefly discuss the effects of fluctuations in the laser pa- 
rameters, of finite system sizes, and interactions between 
the atoms. 



where we have defined the matrix elements 

Tx = J ck w* {x)'w{x — \/4:), 

ry{a) = J dy w*{y)cos{4:aTTy/X)w{y), (8) 

and a = qX/A?:. The values of T^ and Ty{a) as a function 
of the depth of the optical lattice Vq are shown in Fig.O 
Both matrix elements are sufficiently large so that the 
above inequality 51 <C A <C can be fulfilled if we want 
to achieve hopping amplitudes = « of 

the order of kHz. We simplify iJias to find 

ZT \ ^ ( -IX 27:iava t i "U ... \ 

mas = 2^ e „ja„_i^j„ + h.c.j 

-A ^ nation,™. (9) 



D. Total Hamiltonian 

The total Hamiltonian describing the configuration 
shown in Fig. |21 is given by iJ = i/iatt + ^^acc + i^ias 
and we will for simplicity assume = = J . In this 
work we will mainly consider small filling factors of the 
optical lattice ri <C 1 with n the average number of par- 
ticles per lattice site and thus only look at one particle 
effects neglecting the interaction terms Un aird Ux ■ Then 
the Hamiltonian can be written as 

H{a) = (e^'^^^^flt _„a„+i + a,\_„^a„,™+i + h.c.) . 

n,m 

(10) 

This Hamiltonian H{a) is equivalent to the Hamiltonian 
for electrons with charge e moving on a lattice in an ex- 
ternal magnetic field B = 2na/Ae where A = UxCiy 
is the area of one elementary cell. In the remainder of 
this paper we will study the properties of H for neutral 



A. Measurement 

One way of experimentally identifying the number of 
energy bands in the case of rational a is to measure the 
time evolution of the particle density n{y, t) in the lattice. 
We assume the lattice to be loaded from a Bose-Einstein 
condensate Q neglect any interaction between the atoms 
and consider only a single atom. The initial wave func- 
tion of the system is assumed to be 

n , m 

with I vac) the vacuum state and Af a normalization con- 
stant. We then turn on the additional lasers to simulate 
a magnetic field and find the density of particles given 

by 

n{y,t) = («'(t)|at,„a„,„|*(i)), (12) 

to be independent of x. If we choose a = 1/r a peri- 
odic interference pattern emerges in the time evolution 
of l^'(i)). Different paths for hopping around in the lat- 
tice interfere with periodic phase relations (cf. Fig. 0^). 
The periodicity of interference pattern is determined by 
r and repeats itself after exactly r lattice sites as can be 
seen in Fig. ^ for a = 1/6. The periodicity is destroyed 
for values of a which are irrational. An example can be 
seen in Fig.^ where a value of a = l/27r is chosen which 
differs by about 5% from a = 1/6. This little change in a 
is sufficient to considerably alter the density of particles 
in the lattice and to destroy any periodicity. 

B. Parameter fluctuations 

As discussed by Hofstadter (1^ the maximum number 
of energy bands shown in Fig. ^ which can in princi- 
ple be distinguished in an experiment, depends on the 
fluctuations in the parameter a. In our case these fluc- 
tuations are determined by the frequency stability of the 
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a) tj b) 

FIG. 4: Particle density: n{y, t) (in arbitrary units) as a func- 
tion of time t and spatial coordinate y. Light (dark) areas 
indicate a large (small) particle density, a) For a = 1/6 a pe- 
riodic particle density with a period of exactly 6 lattice sites 
emerges, b) For a value of a = l/27r the periodicity in the 
particle density disappears and the visibility of the interfer- 
ence fringes decreases. 



lasers and will thus not be significant. In an experiment 
with neutral atoms the resolution of the energy bands 
will rather be determined by the size of the whole sam- 
ple, i.e. by the number 2L/\ where L is the size of the 
whole sample, and by the spatial resolution in measuring 
interference patterns as described in Sec. IIII Al 



C. Interaction effects 
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FIG. 5: Ground state with interactions: Particle number fluc- 
tuations (Tn,m [a), b)] and superfluid parameter |$n,m| [c), 
d)] as a function of n and m for U — 16J < (7c ~ 5.8 x 4J, 
ujt = 0.06 J and chemical potential /ic = 6 J. The plots a), c) 
[b), d)] show the case a = [a = 1/6]. Light (dark) areas in- 
dicate large (small) values of the functions (black corresponds 
to white to 1). 



For large filling n > 1 the interaction terms in the 
Hamiltonian Eq. [3 especially the onsite interactions be- 
come significant and cannot be neglected. We postpone 
a detailed study of these effects to a further publication 
and only include a graph of how the ground state in an 
optical lattice with a superimposed 2D harmonic trap of 
trapping frequency lot 

V{.,y)^M^{.'+y') (13) 

looks like in the presence of the effective magnetic field. 
We use mean field theory in the form of the Gutzwiller 
ansatz (as described in |2flj), and numerically solve for 
the ground state of the system for finite Un = U . In 
Fig. we plot the particle number fluctuations cr^ „ = 

(("-n,m) - {nn,m)'^)/{n„,Tn) with nn,m = oj, „ja„,m and 
the modulus of the superfluid parameter 4in,m = {<in,m)- 
The effective magnetic field alters these properties of the 
ground state significantly in comparison to the case of 
a = 0. The effective magnetic field leads to a decrease 
of the particle number fluctuations and the superfluid 
density in the center of the trap which is typical for the 
onset of a Mott-insulating phase even for U < Uc with 
Uc the critical interaction strength for the transition to 
the Mott insulator. A similar behavior has already been 
found in |2lj . In addition we find from the numerics that 
interference effects lead to point like increase/decraese in 
1 01 and (T in several of the lattice sites. 



IV. CONCLUSIONS 

In conclusion we have shown that quantum optical 
techniques allow to implement Hamiltonians often used 
to model charged particles subject to an external mag- 
netic field in a lattice. We have shown that all of the 
physically interesting parameter regime can be explored 
by this setup and proposed one method to measure some 
of the most striking features of the fractal energy bands 
of the Hofstadter butterfly. 

The setup we have investigated possesses a lot of pos- 
sibilities for further extensions towards quantum simula- 
tions of strongly correlated systems of charged particles 
like interacting electrons moving on a lattice subject to 
electric and magnetic fields. The major difference and 
at the same time one of the most attractive features of 
the atomic system is a large degree of control that can 
be exerted by quantum optical means. In comparison 
to condensed matter systems the Hamiltonian describing 
the system is very well known and the parameters ap- 
pearing can be controlled and varied over a much wider 
range than is usually the case for strongly correlated sys- 
tems. Also, the time scale over which these parameters 
can be changed is short in comparison to decoherence 
time scales in the system. This allows to study coherent 
dynamical effects that are not easily accessible in most 
condensed matter systems. The detailed investigation of 
such aspects lies beyond the scope of this paper and will 
be dealt with in future publications. 
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APPENDIX A: LASER CONFIGURATION 



sumed to be incident in the xy plane at angles (f)e and 
—4>g, respectively. If a z component of the wave vec- 
tors is avoided the experiment can be done in several 
identically prepared planes of the lattice simultaneously 
which enhances the measurement signal. For a given 
q = ke sin((/)e) + kg sin(0g), A' = (A ujeg)/c ^ kg ~ ke 



fulfilling the requirement that cos(0e) ^ 
we find 



kg cos[ 



We describe the configuration for realizing fii, the re- 
alization of is then straightforward if lasers with suf- 
ficiently different frequencies to avoid interferences be- 
tween r^i and Q2 are used. Two running laser waves with 
Rabi frequencies ^e{g) and corresponding laser frequen- 
cies i^e(g) for driving the transitions |e) ^ \r) {\g) ^ \r)) 
with a large detuning Sr are superimposed. We adiabat- 
ically eliminate the auxiliary internal level |r) and the 
resulting Rabi frequency for the Raman transition be- 
tween |e) and \g) is then given by 



where k. 



(9) 



26 r 

= fce(g){cOs((?!)e(g) ), ± Sin((/)e(g)), 0} 



(Al) 



wave vector of the laser fi, 



(9)- 



the 

Both lasers are as- 



COS0P = 



2qikg - A') ■ 



coa0„ = 



2qkg 



(A2) 



where T = ^{q^ - A'2)(4fc2 - g2 _ 4kg A' + A'2). These 

solutions are only physically meaningful for A' < q < 
^jAkg{kg - A') -I- A'2. Since A' < 1/A the resulting hm- 
itations on possible values of q do not constrain possible 
values of a severely. Only a negligibly small range of 
values of a ~ will not be realizable due to these con- 
straints. We note that if we allow the lasers to have a 
z-component of their wave vector any value of a is pos- 
sible. 
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